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Abstract 



^ \ In this work we re-examine the opacity of the cosmic background radia- 

h : 

tion to the propagation of extremely high energy cosmic rays. We use the 
continuous energy loss approximation to provide spectral modification factors 
for several hypothesized cosmic ray sources. Earlier problems with this ap- 
proximation are resolved including the effects of resonances other than the 
A. 
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I. INTRODUCTION 



With the discovery of the microwave background radiation (MBR) JTJ, it was realized that 
its effect on the propagation of ultrahigh energy cosmic rays (UHECR's) must be considered. 
The first treatments in 1966 by Greisen, Zatsepin and Kuz'min |2||| indicated a sharp cutoff 
for cosmic rays with energies above 5 x 10 19 eV, due to the process 7 +p — > A — > p/nir. 
Physically, the thermal photons are seen highly blue-shifted by the cosmic rays in their rest 
frames so that this reaction becomes possible. 

Recently, several extensive air showers have been clearly observed which imply the arrival 
of cosmic rays of energies above 10 20 eV. In particular, the Akeno Giant Air Shower Array 
(AGASA) experiment recorded an event with energy 1.7 - 2.6 xlO 20 eV J|J|, while the 
Fly's Eye experiment reported the highest energy cosmic ray (CR) event ever detected on 
Earth, with an energy 2.3 - 4.1 xlO 20 eV ||[7|]. Since the GZK cutoff provides an important 
constraint on the proximity of UHECR sources, the origin of such energetic particles has 
became one of the most pressing questions of cosmic ray astrophysics. 

The study of the modification of the spectrum of cosmic rays due to their interactions 
with the MBR en route to us here on Earth, is essential if one is to obtain a deeper un- 
derstanding of the origin of the highest energy cosmic rays. Various treatments of the 
propagation of nucleons in the intergalactic medium exist, including purely analytical ones 



jl3|l , as well as purely numerical ones []1^-|17|. In this paper, we re-analyze the effects of 
meson photoproduction, the dominant energy loss mechanism, using the continuous energy 
loss approximation which is expected to be reasonable for the problem at hand. This allows 
us to discuss, in a general way, the modification of the spectrum of CR's from nearby sources 
due to interactions with the cosmic background radiation. 
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II. PROPAGATION OF COSMIC RAYS IN THE MICROWAVE BACKGROUND 

RADIATION 



The propagation of cosmic rays through the intergalactic medium is described by a kinetic 
transport equation which takes into account three mechanisms for energy losses: redshift, 
pair production, and pion photoproduction. Nevertheless, for nearby sources corresponding 
to propagation times r~3x 10 s years, or distances of less than 100 Mpc, redshift effects 
can be neglected. 

If one assumes that the highest energy cosmic rays are indeed nucleons, the fractional 
energy loss due to interactions with the cosmic background radiation at temperature T and 
redshift z = 0, is determined by the integral of the nucleon loss energy per collision times the 
probability per unit time for a nucleon collision moving through an isotropic gas of photons 
|T8| . This integral can be explicitly written as follows, 

IdE c r», , w v f™ , n(w) 

dw r K(w r ) a(w r ) w r / aw — — (1) 



E dt 2T 2 Jo x ' y ' J Wr /2r w 2 

where w r is the photon energy in the rest frame of the nucleon, and K(w r ) is the average 
fraction of the energy lost by the photon to the nucleon in the laboratory frame (i.e. the 
frame in which the microwave background radiation is at w 3 K). n(w)dw stands for the 
number density of photons with energy between w and dw, following a Planckian distribution 
at temperature T [19|], a(w r ) is the total cross section of interaction, T is the usual Lorentz 
factor of the nucleon, and w m is the maximum energy of the photon in the photon gas. 
Thus, the fractional energy loss is given by 

1 dE _ ckT 
~ E~dt ~ ~2n 2 r 2 (ch) 3 

where k and H are Boltzmann's and Planck's constants respectively, and wq is the threshold 
energy for the reaction in the rest frame of the nucleon. 

The characteristic time for the energy loss due to pair production at E > 10 19 eV is 
t 5 x 10 9 yr |2(| and therefore it does not affect the spectrum of nucleons arriving from 



/*oo 

/ dw r a(w r ) K(w r ) w r ln(l - e -«v/ 2r ^) ( 2 ) 

J WQ 
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nearby sources. Consequently, for nucleons with E > 3 x 10 19 eV (taking into account the 
interaction with the tail of the Planck distribution), meson photoproduction is the dominant 
mechanism for energy loss. Notice that we do not distinguish between neutrons and protons; 
in addition, the inelasticity due to the neutron (3 decay is negligible. 

In order to determine the effect of meson photoproduction on the spectrum of cosmic 
rays, we first examine the kinematics of photon-nucleon interactions. Assuming that reac- 
tions mediated by baryon resonances have spherically symmetric decay angular distributions 
T3| , [T6|| , the average energy loss of the nucleon after n resonant collisions is given by 



1 " / mi, — mi, , 
K(m Ro ) = l--U[l+ \ ( 3 ) 



2 n t=\ \ m R 



i-l 



where denotes the mass of the resonant system of the chain, mu the mass of the 
associated meson, m# Q = ^Js is the total energy of the reaction in the centre of mass, and 



rriR n the mass of the nucleon. It is well established from experiments |2T|] that, at very 
high energies (a/s above ~ 3 GeV), the incident nucleons lose one-half their energy via pion 
photoproduction independent of the number of pions produced. 

In the region dominated by baryon resonances, the cross section is described by a 
sum of Breit-Wigner distributions constructed from the experimental data in the Table 
of Particle Properties 2^] considering the main resonances produced in N'y collisions with 



TrN,7mN, smdirKN final states. For the cross section at high energies we used the fits to 
the high-energy cross section cr to tai(P7) made by the CERN, DESY HERA, and COMPAS 
Groups f22fl . In this energy range, the cr to tai(^7) is to a good approximation identical to 
0tatai(P7)- 

The numerical integration of Eq. (0) is performed taking into account the aforementioned 
resonance decays and the production of multipion final states at higher centre of mass 
energies (y/s ~ 3 GeV). A x 2 fit of the numerical results of equation (Q) with an exponential 
behaviour, v4exp{— B/E}, proposed by Berezinsky and Grigor'eva [|10j for the region of 
resonances, gives: A = (3.66 ± 0.08) x lO^yr" 1 , B = (2.87 ± 0.03) x 10 20 eV with a 
X 2 /dof = 3.9/10. At high energies the fractional energy loss was fitted with a constant, 



C = (2.42±0.03) x 10 8 yr 1 . These results differ from those obtained in | 10| , due to a refined 



expression for the total cross section. From the values determined for the fractional energy 
loss it is straightforward to compute the energy degradation of UHECR's in terms of their 
flight distance. The results shown in Fig. 1, are consistent with those previously obtained 



by Cronin p3], in which just single pion production is considered. Our results therefore 



support the assumption that the energy loss is dominated by single pion photoproduction. 

III. MODIFICATION OF THE COSMIC RAY SPECTRUM 

Let us now turn to the modification that the MBR produces in the UHECR spectrum. 
The evolution of the spectrum is governed by the balance equation 

ON d(b(E)N) nr72sr n 

-m =A ^ + DVN + Q (4) 

that takes into account the conservation of the total number of nucleons in the spectrum. In 
the first term on the right, b(E) is the mean rate at which particles lose energy. The second 
term, the diffusion in the MBR, is found to be extremely small due to the low density of relic 
photons and the fact that the average cosmic magnetic field is less than 10~ 9 G [^| and is 
neglected in the following. The third term corresponds to the particle injection rate into the 
intergalactic medium from some hypothetical source. Since the origin of cosmic rays is still 
unknown, we consider three possible models: 1) the universal hypothesis, which assumes that 
cosmic rays come from no well-defined source, but rather are produced uniformly throughout 
space, 2) single point sources of cosmic rays, and 3) sources of finite size approximating 
clusters of galaxies. In all the cases it has been assumed that nucleons propagate in a 
straight line through the intergalactic medium due to the reasons mentioned above that 
allow us to neglect the diffusion term. 

There exists evidence that the source spectrum of cosmic rays has a power-law depen- 
dence Q{E) = kE"^ (see for instance jnj). With the hypothesis that cosmic rays are 
produced from sources located uniformly in space, with this power law energy dependence 
(which implies a steady state process), the solution of Eq. (§) is found to be 



For the case of a single point source, the solution of equation (f|) reads ]TT| 

1 r c 
b(E) Je 



1 r°° 

N(E,t) = — / Q(E g ,t')dE g (6) 



with 



/ r Ea dE 

t' = t- / — =- 7 

Je b{E) V ; 

and Eg the energy of the nucleon when emitted by the source. The injection spectrum of 
a single source located at to from the observer can be written as Q(E,t) = kE" 1 5(t — t ) 
and for simplicity, we consider the distance as measured from the source, that means to — 0. 
At very high energies, i.e. where b(E) = C E (where C is the constant defined above in the 
discussion of fractional energy loss) the total number of particles at a given distance from 
the source is given by 

or equivalent ly 

N(E,t) w K rv (7 - 1)Ci (9) 

This means that the spectrum is uniformly damped by a factor depending on the proximity 
of the source. 

At low energies, in the region dominated by baryon resonances, the parametrization of 
b(E) does not allow a complete analytical solution. However using the change of variables 

f E 9 dE . . 

t= / 10 
Je b{E) V ; 

with Eg = £(E,t) and di = dE g /b(E g ), we easily obtain, 



oo 



N(E,t) = —j o £(E,t)-i 5(t-t)b[£(E,t)]dt (11) 



and then, the compact form, 



N(E,t) = ^Epb(E g ) (12) 

In our case, where we have assumed an exponential behaviour of the fractional loss energy, 
E g is fixed by the constraint: At — Ei (B/E) + Ei (B/E g ) = 0, Ei being the exponential 
integral [p5|| , and B the constant defined above in the parametrization of Berezinsky and 
Grigor'eva. 

Studies are underway of the case of a nearby source (i.e. within about 100 Mpc) for 
which the fractional energy loss AE/E is small. Writing E g = E + AE, and neglecting 
higher orders in AE, it is possible to obtain an analytical solution and the most relevant 
characteristics of the modified spectrum. Therefore, an expansion of b(E + AE) in powers 
of AE in Eq. (|7]) allows one to obtain, an expression for AE, 



AE < exp 



tb(E) (E + B) 
E~ 2 



ip2 
wna (13) 



To describe the modification of the spectrum, it is convenient to introduce the factor rj, 
the ratio between the modified spectrum and the unmodified one, that results in 

fE + AE\-^b{E + AE) 

Equations ( |13D and (14) describe the spectral modification factor up to energies of ~ 95 
EeV (ps 85 EeV) with a precision of 4% (9%) for a source situated at 50 Mpc (100 Mpc). 
In Fig. 2 we plot the modification factors for the case of nearby sources with power 



law injection (7 = 2.5) to compare with the corresponding results of ref. [ lQfl . The most 
significant features of rj are the bump and the cutoff. It has been noted in the literature 
that the continuous energy loss approximation tends to overestimate the size of the bump 



l4|17|1 . Figure 2 shows how the bump is less pronounced in our treatment. As the bump 



is a consequence of the sharp (exponential) dependence of the free path of the nucleon on 
energy, we attribute these differences to the replacement of a cross section approximated by 
the values at threshold energy flQfl , by a more detailed expression taking into account the 
main baryon resonances that turn out to be important. 
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Another alternative for the particle injection is clusters of galaxies. These are usually 
modelled by a set of point sources with spatially uniform distribution, although in reality 
the distribution of galaxies inside the clusters is somewhat non-uniform. In our treatment 
we shall assume that the concentration of potential sources at the center of the cluster is 
higher than that near the periphery, and we adopt a spatial gaussian distribution. With 
this hypothesis, the particle injection rate into the intergalactic medium is given by 



. (-(T-t 



,2 ' 



Q(E,t) = K -—6(t-T)expl 2 \ dT ( 15 ) 

J -00 V2 7TO- { 2a z j 

A delta function expansion around to, with derivatives denoted by lower case Roman super- 
scripts, 

S(t — T) = 5(t - t ) + 6®(t - t )(T - t ) + ±6<®(t - t )(T - t ) 2 + ... (16) 

leads to a convenient form for the injection spectrum, which is given by, 

Q(E, t) = kE" 1 [ S(t' ~to) + ^ (t - to) + ^ 5^ (t'-t ) + ...) (17) 

From Eqs. @ and (0), it is straightforward to compute an expression for the modification 
factor, 

E~1h(E\ ( rr 2 A 2 p- 2B / E s n A A 4 p~ iB / E 3 } 

" - i 1 + - A ^ F ^ + - Al *— F ^ + °< 6 >} < 18 > 

where 

F 1 (E g ) = 2B 2 E g 2 + (2 - Z^BE; 1 + (1 - 7 ) 2 

and 

F 2 {E g ) = 2AB A E g i + (4 - 50j)B 3 E; 3 + (35 7 2 - 25 7 + 8)B 2 E; 2 
+ (-10 7 3 + 20 7 2 - 15 7 + 4) BE" 1 + (1 - 7 ) 4 

The modification factor for the case of extended sources described by Gaussian distribu- 
tions of widths 2 and 6 Mpc at a distance of 18.3 Mpc is shown in Fig. 3. These may be 

8 



taken as very crude models of the Virgo cluster ||26|| , assuming that there is no other signifi- 
cant energy loss mechanism for cosmic rays traversing parts of the cluster beyond those due 
to interactions with the cosmic background radiation. 

A power law injection (7 = 2.5) was used again, as in the pointlike case. The curves can 
be understood qualitatively. Both peaks are in about the same place, occurring around the 
onset of pion photoproduction. The narrower curve, corresponding to the broader source dis- 
tribution, reflects the losses suffered by the more remote part of the distribution in traversing 
a greater distance to us. We note that the modification factors for nearby extended sources 
can appear similar to those for more distant narrower sources. 

IV. CONCLUSIONS 

We have presented a recalculation of the effects of the cosmic background radiation on 
the propagation of cosmic rays using the continuous energy loss approximation first used 
by Berezinsky [18|] and our best knowledge of 7-nucleon cross sections. This approximation 
is important as it is the only way to allow semi-analytical treatments in a wide range of 
situations. 

While earlier calculations within this approximation suffered from problems described in 
the main text, these do not appear with our improved treatment of cross sections, and the 
results are in reasonable agreement with previous calculations using Monte Carlo methods 
without the continuous energy loss approximation. 

The eventual observation of the GZK cutoff is of fundamental interest in cosmic ray 
physics, providing constraints on the distance to sources of ultrahigh enery cosmic rays. 
The future Pierre Auger Project [p7| should provide a decisive test of the ideas discussed in 



this paper. 



9 



ACKNOWLEDGMENTS 



We are grateful to Huner Fanchiotti and Carlos Garcia Canal for illuminating discussions, 
and to James Cronin for stimulating our initial interest in this subject. We would also like 
to thank Carlos Feinstein for helpful information on the morphology of the Virgo cluster. 
This work was partially supported by CONICET, Argentina. L.A.A. thanks FOMEC for 
financial support. 



10 



REFERENCES 

[1] A. A. Penzias and R. W. Wilson, Ap. J 142, 419 (1965). 
[2] K. Greisen, Phys. Rev. Lett. 16, 748 (1966). 

[3] G. T. Zatsepin and V. A. Kuz'min, Pis'ma Zh. Eksp. Teor. Fiz. 4, 114 (1966). 

[4] S. Yoshida et aL, Astropart. Phys. 3, 105 (1995). 

[5] N. Hashashida et aL, Phys. Rev. Lett. 73, 3491 (1994). 

[6] D. J. Bird et aL, Phys. Rev. Lett. 71, 3401 (1993). 

[7] D. J. Bird et aL, Ap. J. 441, 144 (1995). 

[8] V. S. Berenzinsky and G. T. Zatsepin, Yad. Fiz. 13, 797 (1971). 

[9] F. W. Stecker, Phys. Rev. Lett. 21, 1016 (1968). 
[10] V. S. Berezinsky and S. I. Grigor'eva, Astron. Astrophys. 199, 1 (1988). 
[11] F. W. Stecker, Nature 342, 401 (1989). 

[12] F. A. Aharonian, B. L. Kanevsky and V. V. Vardanian, Astrophys. Space Sci. 167, 93 
(1990). 

[13] J. P. Rachen and P. L. Biermann, Astron. Astrophys. 272, 161 (1993). 
[14] C. T. Hill and D. N. Schramm, Phys. Rev. D31, 564 (1985). 
[15] S. Yoshida and M. Teshima, Prog. Theor. Phys. 89, 833 (1993). 
[16] F. A. Aharonian and J. W. Cronin, Phys. Rev. D50, 1892 (1994). 
[17] S. Lee, Report No. FERMILAB-Pub-96/066-A (unpublished). 
[18] V. S. Berenzinsky, Yad. Fiz. 11, 399 (1970). 

[19] Concerning the density of photons, data from the Far-InfraRed Absolute Spectrometer 

11 



(FIRAS) of the Cosmic Background Explorer (COBE) shows that the spectrum of 
MBR is that of a blackbody of temperature T = 2.73 ± 0.06 K, with no deviation from 
a Planckian spectrum greater than 0.25 % of the peak brightness and a limit on the 
Compton y-parameter of y < 0.0004 (95 % CL). See, R. A. Shafer et at, Bull. Am. 
Phys. Soc 36, 1398 (1991); H. P. Gush, M. Halpern and E. H. Wishnow, Phys. Rev. 
Lett 65, 537 (1990). 

[20] G. Blumenthal, Phys. Rev. D 1, 1596 (1970). 

[21] I. Golyak, Mod. Phys. Lett. A7, 2401 (1992), and references therein. 

[22] Particle Data Group, L. Montanet et al, Phys. Rev. D50, 1173, 1335 (1994). 

[23] J. W. Cronin, in Cosmic Rays Above 10 19 eV, Proceedings of the International Work- 
shop, Paris, France, 1992, edited by J. Cronin, A. A. Watson, and M. Boratav [Nucl. 
Phys. B (Proc. Suppl.) 28B, 213 (1993)]. See also reference p . 



[24] P. P. Kronberg, Rep. Prog. Phys. 57, 325 (1994) 

[25] M. Abramowitz and I. A. Stegun, "Handbook of Mathematical Functions" (Dover, New 
York, 1970). 

[26] B. Binggeli, G. A. Tammann and A. Sandage, Astron. J. 94, 251 (1987). 

[27] The Auger Collaboration, Pierre Auger Project Design Report, 1995, Fermi National 
Accelerator Laboratory (unpublished). 



12 



Figure Captions 



Figure 1. Energy attenuation length of nucleons in the intergalactic medium. Notice 
that, independently of the initial energy of the nucleon, the mean energy values approach 
to 100 EeV after a distance of « 100 Mpc. 

Figure 2. Modification factor of single-source energy espectra for different values of 
propagation distance and power law index 7 = 2.5. 

Figure 3. Modification factor of extended-source energy spectrum for a propagation dis- 
tance ~ 18.3 Mpc and power law index 7 = 2.5. Solid line stands for a gaussian distribution 
of width a = 2 Mpc, while dashed line a width of a = 6 Mpc. 
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